Signature of the electron-electron interaction in the magnetic field dependence of 
nonlinear I-V characteristics in mesoscopic systems. 
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We show that the nonhnear I-V characteristics of mesoscopic samples with metaUic conductivity 
should contain parts which are linear in the magnetic field and quadratic in the electric field. These 
contributions to the current are entirely due to the electron-electron interaction and consequently 
they are proportional to the electron-electron interaction constant. We also note that both the 
amplitude and the sign of the current exhibit random oscillations as a function of temperature. 
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According to the Onsager relation the linear conduc- 
tance G(H) of a conductor measured by the two-probe 
method must be an even function of the magnetic field 

G(H) = G(-H) (1) 

Eq.l is a consequence of general principles: the time re- 
versal symmetry and the positive sign of the entropy pro- 
duction. Therefore it holds in all conductors. It is pos- 
sible however that that nonlinear I-V characteristics of 
conductors contain parts odd in H. In particular, one can 
have contributions to the total current through a sample 
which are linear in H and quadratic in the voltage across 
the sample V. 

I^^i) = aV^H (2) 

Since H is an axial vector and the current is a po- 
lar one, the coefficient a can be non-zero only in 
non-centrosymmetric media. In the case of bulk non- 
centrosymmetric crystals terms in I — V characteristics 
that are linear in H have been investigated both theo- 
retically, the using Boltzmann kinetic equation, and ex- 
perimentally (See for example In the case of chiral 
carbon nanotubes a classical theory of this effect was dis- 
cussed in . 

In this article we study this effect at small temper- 
atures and in mesoscopic disordered samples where all 
possible symmetries are broken. In this situation all elec- 
tron transport effects are of a quantum interference na- 
ture. The theory of nonlinear characteristics of meso- 
scopic metallic samples was developed in the approxima- 
tion of non-interacting electrons ^|^. It is important, 
however, that in this approximation a = and magnetic 
field dependence of the I — V characteristics is an even 
function of H. Therefore the coefficient a in Eq.2 should 
be proportional to electron-electron interaction constant 
/3, which is defined by the interacting part of the electron 
Hamiltonian 
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Here v is the electron density of states. Thus, in princi- 
ple, by measuring the current in Eq.2 one can measure 
the electron-electron interaction constant /3. 

Let us consider a sample of two-dimensional geometry 
shown in the insert of Fig.l and assume that the magnetic 
field is perpendicular to the plane and that the charac- 
teristic size of the sample L ^ Z is much larger then the 
electron elastic mean free path I. At low temperatures 
the main contributions to both mesoscopic fluctuations of 
the conductance SG = G — (G) and the nonlinear current 
Eq.2 are due to electron interference effects. As usual 
in such situations a is random sample specific quantities 
with zero average (a) = 0. To characterize a one has to 
calculate the variance (a^). Here the brackets () denote 
averaging over realizations of a random white noise scat- 
tering potential characterized by a correlation function 
{u{r)u{r)) = j^6{r — r'), where m is the electron mass. 

A diagrammatic verification of Eq.l in the situation 
when interference corrections to conductance are signif- 
icant is not entirely trivial and involves calculation of 
diagrams shown in Fig. la, b J^. We will use a standard 
diagram technique for averaging over random realizations 
of the scattering potential p| . To verify that the correla- 
tion function of conductances (5G(H)(5G(0)) is an even 
function of the magnetic field one has to take into ac- 
count contributions to the correlation function of " Diffu- 
son" and " Cooperon" propagators which are represented 
by the ladder diagrams in Fig.l a,b. (The "Diffuson" 
is the ladder part of the diagram in Fig. la where the 
direction of arrows of the electron Green function is an- 
tiparallel, while the "Cooperon" is the ladder in Fig. lb 
with parallel directions of the arrows. They contain parts 
linear in the magnetic field, which are equal in magnitude 
and of different signs. Thus these contributions cancel. 
To verify the fact that a = in the approximation of 
non-interacting electrons one has to calculate diagrams 
shown in Fig.lc. In this case the linear in H the Diffuson 
and the Cooperon contributions will cancel each other as 
well. This fact is also quite obvious in the framework of 
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Landauer scheme of calculations of the conductance. 
In linear in (3 approximation the variance 
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is given by diagrams shown in Fig. Id. Here F = HD/L^ 
and D jl? is the inverse lifetime of an electron in the sam- 
ple, A is the area of the sample, and D — vpljl is the 
electron diffusion coefficient. Eq. 4 is valid at eV <^ F 
and $ <C $0, where $ = HA and $0 are the magnetic 
flux through the sample area and the flux quanta respec- 
tively. For simplicity we consider the short range e-e 
interaction described by Eq.3. In this case the Hartree 
term is twice larger than the exchange term and we con- 
sider only the diagram shown in Fig. Id. 

In the case of the Coulomb interaction between elec- 
trons and at high electron densities we have 
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where ru is the screening radius ( e^r^/D ^ Ijkpl). 

The existence of the effect described by Eq.4 is con- 
nected to the fact that in the presence of a voltage V 
across the sample there is a part of the local density 



5n{Y) - VH 
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which is proportional to H and V |8j . Existence of fluc- 
tuations of the density of the form Eq.6 is a consequence 
of the fact that linear in H parts of the "Diffusion" and 
"Cooperon" digrams shown in Fig. Id do not cancel. For 
our estimate it is enough to consider only " Diffuson" con- 
tribution. 

We would like to stress that the effect described by 
Eq.4 is quite different from conventional effects in bulk 
crystals which can be described by the Boltzmann ki- 
netic equation 0, Q ■ The leller effects are determined by 
relaxation processes in materials with complicated band 
structures and for this reason they are proportional to 
the the relaxation rates (or proportional to 0^), while 
Eq.2 is proportional to (3. 

The qualitative explanation of Eq.4 is the following. 
The mesoscopic fluctuations of the current density inside 
the sample are due to random interference of electron 
waves traveling along different diffusive paths. Though, 
the total current through the sample should be an even 
function of H the local current densities contain a part 
proportional to H. For example, a part of the current 
density proportional to H can be characterized as a " Hall 
current density" . To avoid confusion we would like to 
mention that this "Hall component" is connected with 
the electric fleld in highly non-local way and has a ran- 
dom direction. By the same token in random system 
there is a component of electron density described by 
Eq.6. We note also that density fluctuations Eq.6 are 
different from Friedel oscillations in disordered samples. 



which are even function of H. In the Hartree approxima- 
tion there is an additional scattering potential 



5ue{v) = —5n{v) 
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associated with the fluctuations of the electron density. 
Thus we can write an expression for a total current in 
the form 



I = G[V,T, H,{6u,{v,H)}]V 
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where the nonlinear conductance G — Gd + 
SG{V, T, H, {Sue{r, H)}), generally speaking, depends on 
the realization of Me(r) via the corresponding dependence 
of the mesocopic part SG. Here Gd = e^Dv is the Drude 
conductance. 

The sensitivity of the sample conductance to a change 
in the scattering potential (5ue(r) has been considered 
in Generally speaking, the mesoscopic part of 

the conductance 5G depends on all spatial harmonics of 
(5Me(r). However, the main contribution to the change 
of the conductance comes from zero harmonics of the 
potential 



5ue{V, H) = y" ((5n(r, V, H) 



Sn{r,V,0))dr (9) 



This can be verified by making calculations similar to 
those in . This is also related to the long range charac- 
ter of the correlation function of the part of the electron 
densities, which are proportional to F For exam- 
ple, in the 2D case the correlation function described by 
digrams shown in Fig.lc has the form 
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Expanding Eq.8 {SI — {dSG / due)ueV) in terms of Ue 
and taking into account that in the main approximation 
u and SG are uncorrelated we get 

{lli))^{{dG/du.f)(ul)V' (11) 

According to [ui {{dG/duf)) = {e^/hf/T^. Calculat- 
ing the correlation function 

^2 ^ ^2 |gy|2 



{S{u,{H)r) = 



(12) 



A^ $0 r 

we arrive at Eq.4. 

At this point we would like to mention that on a qual- 
itative level the effect considered above can be also de- 
scribed in the framework of Landauer-Buttiker scheme. 
To do so one has to combine results of [13 and 0. 

Eq.4 is valid at small temperatures F ^ T. At finite 
temperature the quantity {a'^{T)) decreases with T. At 
T>F 
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We stress that the temperature dependence of a{T) 
is non-monotonic: a{T) exhibits random osciUations in 
magnitude and sign, superimposed on the average decay. 
One can see this by calculating the quantity 

(a(T)a(O)) ^ {a^m^ (14) 

Note that Eqs.13,14 have the same temperature depen- 
dence, which is impossible without oscillations of the sign 
of a(T) Q. 

In the case of high magnetic field $ > $0, (but still 
eV <C r), the part of the current which is the asymmetric 
in H and quadratic in V exhibits random oscillations as a 
function of $. These oscillations, typical for mesoscopic 
systems have characteristic period $0 and the amplitude 

((/(H)-/(-H)f)=/3^(^)^^(^)^ (15) 

Finally, we mention that there are no mechanisms con- 
tributing to Eq.2 other than the mechanism considered 
above. For example, at finite V there is a new channel of 
electron transmission through the sample when an inci- 
dent electron is transmitted into two electrons and a hole. 
The probability of such a process has a component which 
is linear in H. It's magnitude can be estimated in a way 
similar to the estimating the electron- electron scatter- 
ing rate of quasiparticles in a uniform Fermi liquid. As 
a result, it is proportional to V^. Thus the magnitude of 
the asymmetric-in-_ff part of the current associated with 
such process is proportional to V^H. 

Recently discussed above effect was observed experi- 
mentaly in GaAs quantum dots 16]. 

During the preparation of the manuscript we became 
aware of similar unpublished results by D. Sanchez and 
M. Biittiker 
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FIG. 1: Solid lines correspond to electron Green's functions, 
thin dashed lines correspond to the correlation function of the 
random scattering potential {u(r)u(r')) , thick dashed lines 
correspond to the electron-electron interaction l3/i^. Sym- 
bols E, J, n correspond to electric field, current density and 
electron density respectively. Diagrams a) and b) describe 
the correlation function (G(H)G(O)). Parallel and antiparal- 
lel directions of arrows in the ladder parts of those diagrams 
correspond to the Cooperon and Diffusons respectively. The 
diagram c) describes the mesoscopic fluctuations of the part 
of the electron density which is proportional to the electric 
field squared in the approximation of noninteracting electrons. 
The diagram d) describes Eq.lO. The diagram e) describes 
Eqs.4, 13, 14. The insert shows a schematic picture of the 
sample. 



